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ABSTRACT 


The effect of various temperature boundary conditions 
andsaitrerent anlet toroutlet ratios on the load carrying 
capacity of a-plane slider bearing is analysed. The 
Sout LoOUCcTOneOr tne Inertia sterms, and the ettects of 
convection and dissipation are also examined. The lubricant 
is assumed to be incompressible, and the variation of 
viscosity with temperature is taken into account. From the 
governing equations, a set of reduced equations 1S first 
obtained by nondimensionalizing and carrying out an order 
of magnitude analysis. These reduced equations, together 
with the boundary conditions, are then transformed by fee 
o£ the stream function, and the resulting equations are 
solved numerically. Using expressions developed, the load 
Capacity, atag, mass fiow, and average Llow temperature are 


evaluated. 
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CHAPTER I 
STATEMENT OF THE PROBLEM 


The purpose of the present investigation, is to 
study the effects of various temperature boundary conditions, 


on the load-carrying capacity of a plane slider bearing. 
DEFINITION OF PROBLEM 


Up to now, not much work has been carried out on 
methods for fnerese ce the load-carrying capacity of slider 
bearings. Some of the work done by Tahara [1] and a few 
others, have been concerned with removing the heat generated 
in a lubricating film, by use of a coolant forced through 
the bottom of the bearing pad. The reasoning was that 
cooling the lubricant would not only prevent overheating, 
but would also increase the load-carrying capacity of the 
bearing. The fact that cooling the lubricant increases its 
viscosity is well known, but it should also be borne in mind 
that a more viscous lubricant would also tend to increase 
the dissipative effects in the bearing. We can see therefore 
that the main factor contributing to the load capacity is not 
simply the bulk flow temperature. 

It is possible to obtain a better understanding of the 
contributions to load capacity, if we consider the physical 


aspects of the operation of a slider bearing. The main 


sanot2bn30 PN Ee 
swatisaed rei ta ‘eagiq & a ete 
pr 


Go fuo betwies’ need fer, stow oud tort snes oo 
» tebi fe Qos tosses puivxdnocieo. site gotessstont or 
wet ep bas Lt) eerie? Yo Sook teew edd Qe eeipe) 
bertsishSp Basi oid palvemer tetv HEALED Abed ¥en | 
dpinerds -beotc? ganscas. 4 Fo sauce ries! pmisheridet a us 
gait esvcoeineeks: alt. beg paiteed edd ideal at av 
<preiseonieve Jrevety yinn, ton: Biluow, see pbtdecl sity susttel’e 
sna to vajuscess PiiVssesrSbol. sds caseiend dels’ a 
asi eomaarans dhsotadul aid ehbideap seme 36% _ 
bellu mt oiryod o¢. cele Fuori a2 Sud vail ei Bi ¢ 
Sef oni oF tues Oele bloow jnshsauas “guobety a & pe 
Stohemses “Se uso Sw -ynlased edt as esuatte evideginere e “ 
dom Bk yeiveaes Seok aris of pa taedi aioe edoes ie 
prscore i 
cae 


ant Fo Sonera tees iadted 6 pacha lee 


_ 


nequirement wor the operation of any bearing, 1s that the 
bearingemiisuipe able ‘to support a, load. In the case of 
fluid film bearings, this requirement means that the 
lubricant must be under pressure. For shafts operating 
at very low speeds, and under heavy loads, this pressure 
is usually supplied from outside; and these bearings are 
known as externally pressurized or hydrostatic bearings. 
In the case of bearings operating at higher speeds, and 
under average loads, the pressure is usually supplied by 
the revolving motion of the shaft which drags lubricant 
into the bearing. These types of bearings are usually 
referred to as hydrodynamic bearings, and are the ones 
Cliewicerest in this sseudy . 

When the lubricant enters the bearing, the tapered 
shape of the pad surface (See Fig. 1) forces the fluid close 
ROmLGuLOnfOW parable l) cor the nclinedusun trace. nln reaction 
to this constraint, the fluid exerts a pressure on the bearing 
pad. This pressure builds up to a maximum around the 
middle of the bearing then decreases to the ambient pressure 
atetie bearing OuGret. 

The amount of load a bearing can carry is influenced 
by the inlet to outlet ratio, h./ho- It has been shown [2] 
that the load capacity rises to a maximum for h,/he = a2 2 


then falls gradually. For h,/h, fixed, the load which the 
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bearing can support is determined by the thickness of the 
lubricant film. For the same load the thicker the film 

the larger the load the bearing can carry. This: is because 
it will require a larger load on the bearing in order for 
Pheer ill ch ickiess lO retaai ort s Ori ginalevalue, 

OUr dattempes at Mncreasing the load capacity should theretore 
be concerned with increasing the thickness of the lubri- 

cant film for a given load. This can be accomplished 

MOSG ecliectively by cooling the siider, or iby, cooling the 


bearing pad. 


INERTIA EFFECTS AND TEMPERATURE VARIATION 


IN A SLIDER BEARING 


The Reynolds equation is based on the assumption that 
the inertia forces are negligible compared to the viscous 
forces [3]. The importance of the inertia terms relative to 
the viscous terms in the equation of motion, can be charac- 
terized by a dimensionless parameter referred to as the 
modified Reynolds number, Few 

In the operation of most bearings, ae is very small, 
therefore the inertia effects could be easily ignored. 
However, for bearings operating at high speeds, or where the 
Viecosttysotethe Lubricant is very lowe (erg. Lor gas bearings), 
Re~ may reach values. near to or exceeding unity. Thus, in 


such a case, the inertia forces become comparable in magnitude 


to the viscous forces and may no longer be ignored. If the 
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5 
Reynolds number becomes greater than a certain critical value, 
turbulence will develop, and the governing equations will 
have to be correspondingly modified. 

Another assumption made in classical hydrodynamic 
Jubrication, is that the antares pigctatlel eaey@l aligl. Saleks 
lubricant film is negligible... The viscosity is. then taken 
as being constant throughout the film, and equal to the 
InleLeolleviscosity.  Thiceresults in theyuncoupling of the 
momentum and the energy equations, which greatly simplifies 
the solution of the problem. Under most operating conditions 
though, the energy generated by viscous dissipation causes 
a significant rise in the lubricant film temperature. This 
results in a corresponding decrease in viscosity. Hence, 
the solution assuming viscosity at the inlet oil temperature 
will be in error. Wilcock [4] and Rosenblatt [5] conclude 
from their experiments, that in fact the average value of 
viscosity corresponding to.that of the outlet oil should be 
taken. Several authors have attempted to include the 
VISCOSity Variation, by allowing the viscosity to, vany lin 
the X direction only. This results in a much simpler form 
of the momentum and the energy equations. 

The temperature distribution within the bearing, 
depends largely on the thermal boundary conditions imposed 
on the film. If both boundaries were taken as adiabatic, all 
the heat generated by viscous dissipation would have to be 
Carried out by the fluid. The effect of this would pe to 


decrease the average lubricant viscosity, thus reducing the 
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load carrying capacity of the bearing. If the boundaries 
were not adiabatic though, some of the heat generated would 
be transferred across the bearing faces, and the resulting 


decrease in load carrying “capacity would not be as great. 
REVIEW OF RELEVANT LITERATURE 


The equations of motion for the slider bearing, yield 
two non-linear partial differential equations; and over the 
years, several authors have investigated the influence of 
the non-linear inertia terms using approximate methods. 

Slezkin and Targ [6] averaged the inertia terms across 
the film thickness. This allowed the resulting equation to 
be readily integrated, since the inertia terms became a 
function of x onty.- Other ‘authors [7,8] have employed this 
technique to obtain solutions for various bearing 
configurations. 

Kahlert [9] assumed that the inertia forces were 
small compared to the viscous forces. He then included the 
influence of the inertia terms by a small correction to the 
results obtained from the purely viscous consideration. 

Snyder [10] obtained a more exact solution, by taking 
interaccount tthe variation tof “the anertirateifects “across gehe 
film kasmwe ll eaAsiuinethesdirvection of tilow ws? elm this. methods, the 
stream function is expressed as a power series in 6, which 
ieWarfunction ‘of film thickness; while ‘thercoecificients@are 
assumed to be functions of 1, which in turn is a function of 


both x and y. These coefficients are obtained by solving a 
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set of differential equations. In his solution, Snyder solved 
for only the first two terms of the power series. More 
recently, Rodkiewicz and Anwar [11] found that it was necessary 
LOuretain at least four terme of the series expaneion, in 

order to improve the accuracy at higher modified Reynolds 
numbers. 

In all of the previous references, the solutions were 
obtained assuming constant viscosity. Other researchers have 
EAkenetncoraccount thetvariation of “the Jubricant viscosity, 
for various boundary conditions on the film temperature. 
Hunter and Zienkiewicz [12] considered two different boundary 
conditions on the temperature. 

(1) Both boundaries having the same temperature. 

(2) Both boundaries adiabatic. 

In their analysis, they neglected the inertia terms, and also 
the lateral convection term in the energy equation. 

Synder [13] studied the overall temperature variations 
across the lubricant film. The slider surface was taken as 
an isothermal surface, and the stationary surface as adiabatic. 
Feom Nis analysis, he “concluded thatthe consequences of the 
Variations of viscosity cannot be neglected for high Prandtl 
number lubricants. The most complete investigation so far 
has been carried out by Hahn and Kettleborough [14], who 
used matrix methods. They retained the inertia terms in the 
momentum equation, and the convective terms in the energy 
equations) Also, the viscosity of the lubricant was@regarded 


as a function of temperature and pressure, and the density a 
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function of temperature. The temperature boundary conditions 
were obtained by equating the heat fluxes at the solid-liquid 
interfaces. They concluded that it was desirable to cool the 
bearing surfaces, since this would increase the load carrying 
Capacity of the bearings. It was also found that the 
temperature at the moving surface was reasonably constant, 


even though the heat flux varied considerably. 
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CHAP DE helt 
THE GOVERNING EQUATIONS 
GENERAL FORM OF THE EQUATIONS 


The equations governing the steady, laminar flow of 
an incompressible lubricant sthrotighuan inmtinite slider 
bearing, as shown in Fig. I, are the momentum equations in 
Roce Spahtimy directions, (Has (2.1) sand (252) respectively. 


thesenergy equation, Eq. (2.3), and the continuity equation, 


Eq. (2.4). In order to solve these equations, it is also 
necessary to specify the boundary conditions as presented 
Asta CWS eae | Cun) be 

For our problem, the boundary, conditions om velocity 
are obtained by satisfying the conditions of no-slip at the 
pad and at the slider; while the boundary conditions on the 
temperature, are obtained from the compatibility conditions 
at the said two interfaces. 


The general form of the equations may now be listed 


as follows: 


Momentum equations 
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Continuity equation 

du OV 

ax ages 0 (2.4) 
Boundary conditions: 

y=0; u=U, v=o 

y=h; u=0, v=o 

20%) 
y=0; T=T. 
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The equation used for the film thickness is 
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ALaL 
REDUCTION OF THE GENERAL FORM OF THE EQUATIONS 


BY AN ORDER OF MAGNITUDE ANALYSIS 


Because of the large number of terms in the Governing 
Potatvonce (| 2. be chroughi (2.4), 81 tewouldbesvery dtrriculteco 
Obtalnma satistactory Solution. In orderlto simelatyucne 
problem, the given equations will be nondimensionalized, and 
an order of magnitude analysis carried out. It would then 
be possible to retain only the most important terms, and 
obtain solutions for the derived set of reduced equations. 

Before proceeding with the nondimensionalization, a 
list of some of the usual assumptions made is in order: 

(Lj shhesubricanc Gilme@isavery thin ise. h/, eels 

(2) The density, thermal conductivity, and specific 
heat of the lubricant are constant. 

(3) The effects of thermal and elastic distortions 
are negligible. 

(4) Cavitation effects may be ignored. 

(5) The temperature distribution at the inlet is 
linear. 

Selecting appropriate reference quantities, the 


variables can be expressed as follows: 
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By substituting these transformed variables into 


the governing equations, and with some rearranging, it is 


possible to compare the relative magnitudes of the various 


terms. 
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transforming the x derivative, we obtain 


Oe 
ay 


a] ar fe 


Sl 


lie 


6), 


h 


et (h,-ho) ie >) By Alls 


— (is 
( af 


= ho) /i 


~(h,/h, ~L)h o/h 


oO : 
ay fintriatones’ aneig yn 


oe ae, | 


ai #1) .onbenak=ties shee aac bas ve 


nape yalnse 


_ Y 


bnokesy ofg Yo eabustaven svidsteos ous STH ip 
° ~os 


_ 


- 


a pA ietse isfazs coaU 


i 


Similarly for the y derivative, 
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For the continuity equation, we get 
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The momentum equations become 
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and the energy equation yields 
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One of our earlier assumptions was that the 
Pipi icatingwer iim Was’ Very thin, 16.0) b/im << ee Prom iis 
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these facts in mind, it is possible to retain only those terms 
with the larger coefficients. We are thus able to obtain a 
set of reduced equations, which are to be solved for a given 


set of boundary conditions: 


For the momentum equation in the y direction we get 
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Now both Ea OSs. u ake and h are of order [1]. 
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Therefore 2B is of order Re (—) (2s) 


which is very small, even for modified Reynolds numbers as 


large as 1. We see then that for our problem, the pressure 
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gradient in the y direction is not important, and can be 


tCakenyto equal 0. 


Thus p Ls not va, function of Y;, but of x alone. 


equations when reduced give 
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In the above, the dimensionless parameters are: 
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and the boundary conditions become: 
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GamAS?. eAxateba0 (2.20) 
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The expression derived for the dimensionless film 


thickness is, 
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TRANSFORMATION OF THE GOVERNING EQUATIONS 


In order to reduce the number of dependent variables, 
thus simplifying computations, we can resort to using the 


stream function, Writing, 
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The continuity equation is identically satisfied, and ae 


represents the mass flow through the bearing. 
Expressions (2.15) and (2.16), (see Appendix A), and 


Equations (2.19) now become 
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Expressions (2.20) through (2.22) however, remain the 
same. 

It can be seen from the above equations and the 
boundary conditions, that in general, the solution to our 


problem depends on the following dimensionless groups: 
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CHAPTER 211 
SOLUTION OF THE TRANSFORMED EQUATIONS 
INTRODUCTION 


Because of the non-linear nature of the momentum 
Equation (2.24), and the coupling with the energy Equation 
(2.25) through the viscosity, the solution of the problem 
can present some difficulty. For modified Reynolds numbers 
much less than unity [11], solutions were obtained neglecting 
the inertia terms; while for modified Reynolds numbers in 
the vicinity of unity, solutions were obtained including the 


inertia terms. 


SOLUTION OF THE MOMENTUM EQUATION NEGLECTING 


INERTIA TERMS AND FOR CONSTANT VISCOSITY 


Neglecting the inertia terms, and for constant 


viscosity, the momentum Equation (2.24) becomes 


935 = 
= = ie ele (S215 
dy Y falbie 


Upon integrating the above equation three times with respect 


to y, we obtain the following expression 
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determined from the following three boundary conditions. 
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allows us to obtain an expression for the pressure gradient, 
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where D, (x) and D, (x) are given by the expressions of 


Bavationy (3.6) below 
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The dimensionless mass flow, Gee is found by 
substituting the expression for dp Of hQitar1 On 35.5))e, 


dx 


IncvOernesintegralpErquationy (220), yielding 


oS we D(z) dk / po D(X) dx (3.7) 
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SOLUTION OF THE MOMENTUM EQUATION NEGLECTING 


INERTIA TERMS AND FOR VARIABLE VISCOSITY 


When the inertia terms are neglected, but variable 


viscosity is preserved, the momentum Equation (2.24) yields. 
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The above equation is integrated three times with 


respect to y, whence the following equation is obtained. 
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The coeffs. Cy (x), Gh Ax) and C(x) of Equation 
(3.9) are determined from the boundary conditions as listed 
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The fourth boundary condition allows us to determine 


an expression for dp in the same form as Equation (3.5). 
dx 


D, (x), and D, (x) now become: 
— _2 —= —_ 
D, (x) =-=}.0/h [A, (x,1) - B. (x) ] 
Ciel) 
te = ee) = _ 2 
Do (x) = By (x)/ h [A, (x, 1) B, (x) ] 
An expression for ve is likewise obtained, and 
retains the same form as it does in Equation (3:7). In 
the above, 
B. (x) = A, (x,1) TA, (x,1) = A, (x,1) ]/1A, (x1) Fs A, (x,1)] 
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and 


Al v7 = be 
Are, ly (= fe tae (ae : dy) dy 
De Pee se 
hy (iy WD ee = 
(3548) 

= A 7 Vaal = = 
A, (x, 1) = es Is 7 dy) dy 

= i ae os 


We are now in a position to evaluate Cy lx) and aCa(x) iy 


2 


provided the viscosity distribution is known. 


SOLUTION OF THE MOMENTUM EQUATION INCLUDING INERTIA 


TERMS AND FOR VARIABLE VISCOSITY 


The momentum Equation (2.24) can be rewritten as 


follows: 
2- -3 - 
J Yee Ee ee aay ea) Se (om) 
a ¥, ak 
om, y eS 
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Integrating Eq. (3.14) we obtain 
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C3 (x) = Q 
(edt lus) 
The pressure gradient, ay has the same form as 
ale 
expressed by Equation (3.5), where now 
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In the above, 


SIGS aval: A,(x,1) [1+ A, (x,1) - Ne ay 
[A,(x,1) - A, (x,1)] 
Bo (x) = 0A, (x,1). [A)(x,1) = A, (x,1)1/1A,(%,1) - A, (x,1)] 
eto)) 
and, 
meceAind Sauron «¥ Cup Meng yerdgydy dy 
Deane fo) OF ane yiey Y 
(3.20) 
FENiconY meen. (2s, ARC eet 
6 u fe) 7 e) Ms 


Thus, providing us with the viscosity distribution, 


would allow C, (x) and oe 2) to be evaluated. 


SOLUTION OF THE ENERGY EQUATION 


For the case of variable viscosity, it is necessary 
to obtain the viscosity distribution, before the stream 
function can be computed. This is accomplished by first 
solving the Energy Equation (2.25) together with the boundary 
Conch Ons 2.2. Jenlsing the, Crank-Nicholson.method ea the 
viscosity distribution is then obtained by interpolation 
from a set of data which had been previously read in. 

In the numerical process, (See Appendix B), the 


derivatives are replaced by the central finite difference 
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approximations, and the terms are rearranged to obtain the 


following equation. 
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where ais! paar ain and th ue are defined - in Appendix B. 


Te sould be noted-that the coerricients a. 


ee b 


a euaKel 
uy i 


ese are known from the solution of the momentum equation 
alone. However, sa 2 in addition depends on the temperature 
at the preceding section. The temperature distribution at 
the inlet is given by the last of Equations (2.21). Using 
these initial values, and marching from section to section 
along the bearing, one obtains the complete temperature 
distribution. 

For the numerical procedure, the square grid presented 


below was employed, with the spacing interval in the x and 


the y directions taken as 0.1 and 0.05 respectively. 


For 2 =.1, Equation (3.21))becomes, 
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Using values of j from 2 to N-l1, the above equation 


yields a linear system of N-2 equations with N-2 unknowns. 
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Known OFF lea] i456 ie ve sk, thei system) cba nquamonse(se23) 
Gan be rewritten in matrix form. “The matrix is 
Chto lagOnalwandsOoleOorder N=Z uel t was SOlvedudtbeculy 7 oy 
using the method of triangular decomposition [15]. After 
T. \=—2 5,20 N-1) was evaluated,= tne temperacure 


distribution throughout the grid was obtained by repeating 


the above procedure for each successive section. 
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METHOD OF PSOLUTION = "STEP eBY< STEP 


In the general case, the solution of our problem 
requiressthatswe solve Eqs.) (2.24) "and) (2.25), stogether with 
ches boundary CONG PE1ONSHases bacelm ins hOs mil 22 1) mratnidan(202,6)) 0. 
Since Eqs. (2.24) and (2.25) are coupled through the viscosity, 
an iterative method is used, which involves a back and forth 
movement between these two equations, until the temperature 
distribution converges to within a prescribed error bound. 
Solutions are obtained for the case in which inertia terms 
are not important, and also for the case when these terms are 
important. For the latter case, the solution neglecting the 
inertia terms is used as a first approximation, then these 
terms are included by use of a correcting procedure. The 


method of solution is as outlined in the following steps. 


(1) Evaluate By ts) and Dy (x) by using Bas. (376). 

(2)> Bind the: mass “flow by use of Eq. (3.7). 

(3) Use Eq. (3.5) to determine the pressure gradient, 
ap 
ax 

(4) Determine the value of the coeffs. Cy (x) and 
C(x) bywusingeeqs.w1(3.4)). 

(5) From Eq. (3.2), evaluate the stream function 
duster bution. 


(6) Obtain the temperature distribution by solving 


Bow (2225) €or, constant viscosity. 
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29 
Evaluate A, (x,1) through A, (x,1) by use of 
Bastian 3s ) 
Use Eqs. (3:12) to find B, (x) and B, (x). 
Determine D, (x) and D, (x) bYFUsSe OLfSEqs< 9 (3411). 


Pind the®mass fllow by °usevof Eq?* (897); 


Use Eq. (3.5) to compute the pressure gradient, 
dp 

dx 

Evaluate the coeffs. (e (x) and C, (x) by use 


OL EGS., (3.40) 

Obtain the stream function distribution by 

BIS eee Vkepr eu eth ehh A 

Determine the temperature distribution by 
SOLVIT EC .s8 (2 20 jee 

If, the convergence criteria on the temperature 


is satisfied after the first iteration go to 


ie ethe "Convergence Coiceria On egtne, cermperature 
1S NOt. Satisiied, atrter the first; teeration. go 
£0 step (7). 


Evaluate the bearing characteristics and print 


Tf inertia terms are not important, “goto 


Te anertia terms are amportane continue. 
Use Ha. .(3.15) to £ind the *correctivon stor the 
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presence of inertia terms. 
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(21) Evaluate A, (x,1) through A 


(Geran ish) aewelet (Savoy) = 


¢ (XL) byewuse sol Eqs. 

(22)5Use"rqs. (8.19) to find B, (x) and B, (x). 

(23) Compute D, (x) and D, (x) byausing eqs ..0(3.15)5 

(24) Determine the mass flow by use of Eq. (3.7). 

(25) Use Eq. (3.5) to find the pressure gradient, 
hence determine the pressure distribution. 


(26 aEveliiatemthe ccs. mmeCu Go mandsC@ aKa) by use 


iM Y 


Of rEqQS. "(SR 1a/)e 

(27) “Obtain thesstreamifunction distribution by “use 
OLreh gras 6) 

(28) If the convergence criteria on the pressure is 
satisfied go to step (30). 

(29) EGOStOSet enim): 

(30yMSClve Bq. (2.26) for the temperature distribution. 

(31) Df Sehetconvergence Ycriteria son’ thestemperature 
LSeno besatistied Vattenethe thirst. ternation, go 
to «step ~(20)s: 

(32) Evaltate the bearing characteristics and print 
results. 


(33) Solution completed. 
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CHAPTER IV 
RESULTS AND CONCLUSIONS 
INTRODUCTION 


In order to understand more fully the various 
factors which contribute to the load capacity, the 
Momentum Eq. (2.24) together with the Energy Eq. (2.25) 
are solved for several different cases, which are as follows: 

(a) The momentum equation neglecting inertia terms, 
together with the energy equation are solved for various 
boundary conditions on temperature. Sample results are 
listed in Tables la through ld, and graphs are presented in 
Pigs 14) cncougn L5. 

(b) Assuming constant viscosity, the momentum 
equation is first solved for various values of modified 
Reynolds number. Afterwards, it is solved together with the 
energy equation, for three different temperature boundary 
conditions, but this time retaining the viscosity variation 
with temperature. Tables 2a through 2c list some sample 
results, and graphs are’ presented in Figs. lo“through 23. 

(c) Considering variable viscosity, the momentum 
equation neglecting inertia terms, together with the energy 
equation, are solved for Various Values of “inlet Horoutler 
retio, hy, Sample results are listed in tables 3a through 


3c, and graphs are presented in Figs. 24 through 35. 
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(d) Neglecting the inertia terms, the momentum 
equation together with the energy equation are solved for 
three different ey aks ratios, and tor? constantere :? wit lesat 
thes same’ time retaining the viscosity variation. with 
temperature. In this case, the results obtained enable us 
to examine the effects of dissipation on the load capacity 
ofethe bearing .t)-Graphs are plottedsins Figsss 36"throughs 46, 
and? sample results are disted in tables 4a through 4d. 

(e) In order to examine the effects of convection 
Cermsmoni they load tcapacity of therbearing). ajmethod cimigdar 
to that in (d) was undertaken, but this time keeping PE 
constant instead of Pe. Graphs illustrating results are then 
plotted in Figs. 47 through 57, while sample results are 


presented in tables 5a through 5e. 


EPEEBCTS OF VARLOUS TEMPERATURE 


BOUNDARY CONDITIONS - NO INERTIA 


ha CiShecartnnoughiiGrp ormithesy velocity idisert boutons or 
Eney £foLlowingmboundary conditions: 7a 1) T = 200°F, A = 1OO°F, 
(a 78) T a 100°F, a =) L007 Feandi Gs 2am) TS 5:0 oa male b= “S00 Re 
All the dimensionless groups have the same values, and these 
boundary conditions were chosen to find out what happens to 
themloadseapacity when the slider is cooled. It “is possible 
to get some idea of what the pressure distributions for ‘the 
three different boundary conditions are like, by examining 


their velocity and temperature profiles. 
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If we first examine the shape of the velocity 
dietrubaivion atevanilous isecti onsjivas Gldiustrated: tmuanylor 
EROQStEs po, ODLOpewe sseemthatrthe profile changes from concave 
inwards at the entrance to convex outwards at the exit. This 
is readily explained, because the pressure gradient, which 
is positive at the inlet, passes through zero. around the 
mid-section of the bearing, then becomes negative at the 
exit. Also, by examining only the inlet and outlet velocity 
profiles of Figs. 4 through 6, we can tell which pressure 
distributions are higher. It is reasonable to expect that 
large positive pressure gradients at the bearing inlet would 
tend to cause backflow, while large negative pressure gradients 
at the exit would cause the velocity profile at this section 
to be more convex outwards. 

These two effects can be clearly seen in Fig. 6. Note 
that the outlet velocity profile becomes more convex outwards 
as the slider is cooled, indicating an increased mass flow 
rate. The temperature distribution for the case T. = 200°F, 
Tae LOOGE wappeare sin Fig. 07, whencesthe cise any temperature 
as the lubricant flows through the bearing is readily seen. 
Plots of the,ptessure distribution for the three cases appear 
in Fig. 8, and agree with what we were led to expect from 
observations just made. 

Fig. 9 presents plots of the shear stress distribution 
at the slider for the three cases being studied. Since T = 

au , then the shear stress at the slider depends on both 
Jy | 
the viscosity and the velocity gradient at the said interface. 
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Examination of Figs. 4 through 6 show that the magnitude of 
the belocity gradient at the slider decreases as the slider 
temperature decreases; but since the viscosity rises quite 
rapidly with decreasing temperature, it is very likely that 
chercolder*stiader®willcausemlarger’ shear stresses. to be 
developed. This is indeed the case as can be seen from the 
euLvVesseplotted Gn Figt, 9 

Note that for the colder slider temperatures, the 
velocity gradient passes through zero, then becomes positive 
towards the bearing outlet. This explains why the shear 
SeLesco Curves Cross, cach Other, and gapproach Zero at) polncs 
at or just  before™the* outlet. 

Figs. 10 through 12 present curves of load capacity 
versus ell Ss for the following cases: (1) oes fixed at’ VUGCr, 
Sader cooltcu. a1) T fixed™at  L00°cr, pad. leaced an 
fixed at 350°F, slider cooled and at the same time pad heated. 
Comparison of these curves show that providing Pe, PE, and 
i are the same, the toad» capacity for each of "the three cases 
Just MencLonea, and Lor the same Beg Tg ratio are quite close, 
although the actual boundary’ conditions differ greatly.” It 
can also be seen that the divergence becomes smaller in the 
range. Uv. = Bal So 2.0, which includes the normaly operating 
range of most bearings. 

Figs. 13 and 14 present curves of drag versus Bags 
POLetieme loot, tWwO Cases listed in the previous paragraph. 


Applying a similar reasoning, it is seen that the difference 
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actual divergence jinvload.capacity} drag) ‘etc. for the ‘same 
fev ratio is due to the differences in viscosity distribution 
which follows because of variations in the temperature boundary 
conditions. Fig. 15 shows the influence of ‘the actual 
temperature boundary conditions on load capacity, for two 
different vas tatkvos.: Tigussscenstiavebotm Varger fand 
smaller values of Bey as tend to cause greater divergence in 
the load capacity values, but this divergence is much more 


pronounced for the smaller oy as Lavior 
INFLUENCE OF INERTIA TERMS 


Figs. 16 and 17 present the velocity distributions 
for the cases where inertia terms are not included and where 
they are included respectively. It is seen that when the 
inertia terms are included, the velocity gradient at the 
Dil eiawus, close, tho. zero Mearttorthe ypadis bute aseti ses igh bly 
greater than zero at a similar point for the inertia-less 
case. Thus the pressure gradient at the inlet when inertia 
terms are included, should be slightly larger than for the 
case of no inertia. However, the difference in the outlet 
velocity profiles is not easily discernible, so we would 
expect the pressure gradients at this section to be quite 
close. 

Temperature distributions for the two cases being 
discussedsare presented) in’ Figs.0 18 andelo respectively. = The 
temperature profiles are quite close, and there are no 


significant differences between the two. It is possible 
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though for both cases, to observe the heating of the lubricant 
as it flows through the bearing. 

Fig. 20 presents curves of the pressure distribution; 
whence it is seen that inclusion of the. inertia terms gives 
a slightly higher pressure Anise tac which in turn yields 
eg argertload-Ccapacity, » Fag m2 presents the. shear Stress 
aise: bat ony for Re = 0.5 and for the cases with and without 
Vhewehertaa terms.) »ltyshould. bes noted. that) although 
inclusion of the inertia terms gives a higher pressure 
distribution, it also causes slightly larger shear stresses 
at the slider. 

Fig. 22 compares the results of the present work with 
those of Rodkiewicz and Anwar [11], and Woodhead and 
Kettleborough [15]. It is of interest to note that the curve 
for the present work is almost parallel to the one obtained 
ie tebieaetlaking, the, resultsnot [(L5],assa Standard,,ut is 
seen that the present results are quite good for small and 
for large values of pel and the error rises to a maximum 
OF) el) stor coe around: 026% 

Pi Gemeticy DGesentsycumvescsshowi ngmthesnt bience OL rhe 
iner tilassteonmswon load Capacity mon dLEpernents temperature 
bouncdany condi tions, while table: 2c) lists ithesactual values 
of the percentage increase in load capacity. It should be 
Noted thatethe largest errors made. in neglecting the inertia 
terms occur when the bearing operates under normal conditions, 
ie. when a SiliO. meron COnctan fj avscocley mand ton ee 
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wiich isa bit less than the value of 8.9% given in [11)% 


EFFECTS OF VARYING INLET TO 


OUTLET RATIO - NO INERTIA 


Figs. 24 through 26 present the velocity distributions 
for the following cases respectively: (i) Bt SS (de) 
hy Sues, ge oee( ic) hy =" 2/55’ Notice’that’ as in increases, the 
backflow at the inlet of the bearing increases continuously, 
while the velocity profile at the outlet becomes more convex 
outwards. These observations would lead us to expect a rise 
in the pressure distribution as hy increases from 1.75 to 
2./5. However, we must first take a look at the temperature 
distributions. 

Figs. 27 through 29 present plots of the temperature 
idstributions for the three values of hy under consideration. 
The respective profiles are quite similar, so that the 
temperature distribution does not seem to have contributed 
much to the differences observed in the respective velocity 
distributions. 

As shown by the curves of Fig. 30, the expectation 
of a continuous rise in the pressure profiles with he Ss OnLy 
Parclys correct.9° What can be"readily seen"is thal for positive 
gradients, the smaller the value of Ae the larger the pressure, 
while for negative gradients the reverse is now true. 

Figs Sl@presents plotsvof the sheer stress disetri— 
bution for the values of he under consideration. It should 


be noted that the magnitude of the shear stress decreases as 
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hy increases. In dimensionless form, the shear stress at 


the slider can be represented as follows: i7 = (U oes : 
~ dy y=0 
Since the slider temperatures are all the same, and the 


sectional gradients ae) for the three values of hy are 


note toorditferents,® then ee eke in the shear stress 
profiles will be influenced mainly by the variations in l/h. 
Since h is Proportronaleto hes then» 1/h.a L/ tes So mhag 
increasing hy should cause numerically smaller values of 
shear stress at the slider. This is.indeed the result, as 
San begseenyfrom the curves of Fig. 31s 

Figs. 32 through 34 plot curves of load capacity 
versus hs for constant Pe, but for different temperature 
boundary conditions; while Fig. 35 compares curves with 
different temperature boundary conditions, but for constant 
Pe. For all these figures, it should be noted that the load 


capacity always rises to a maximum for 2.00 < he Ke 2254 


then falls gradually. 
DISSIPATION BRFECTS — NO INERTIA 


Pigs.as6 throughes8t ploteathesvedocityhdistrm butions 
Lop re ils 275ebn Oseands 20n SSurespectivel yiqwhitesortner 
dimensionless groups have the same values. Upon examining 
these figures, we find that the velocity profiles become more 
convex outwards as PE increases, indicating a rise in mass 
flow. This trend might lead one to expect an increase in 
load capacity also, but this conclusion cannot be made 


without first examining the temperature distributions. 
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If we now examine the temperature distributions 
of Pigs. 939 through 41) wesfind*’ that theidifferences among 
the three are quite pronounced, with the temperature profiles 
rising continuously as PE increases. The pressure gradient 
is directlyeproportional to the» lubricant viscosity) which 
decreases with increasing temperature. Thus the increasing 
mass flow as PE increases is counterbalanced by the drop in 
lubricant viscosity; and whether the load capacity rises or 
falls with PE, will depend on the combined influence of these 
Ewortactors. 

The curves of Fig. 42 show that the pressure profile 
gets smaller as PE increases. So that the reduction in 
viscosity seems to have had a much greater effect on the load 


Capacity, than the small rise in mass flow of the lubricant. 


The shear stress is given by To= i: and since the 
. h 
lubriicantlin contact with thesslider has ajyconstant viscosity, 


then the shear stress at the slider is proportional to the 
velocity gradient at. the said interface. Examination of Figs. 
86 ehrough (38 show thatethesvelocityeqradventeatithe sisder 
are all. negative, and that the magnitude of the gradient at 
any section decreases as PE increases; so that the effect of 
inereasingeausstipationsis to decrease the drageonsthesstider: 
Results for the shear stress are plotted in Fig. 43, and 
agree quite well with the observations just made. 

Figs. 44 through 46 plot curves of load capacity 
versus PE for constant Por but for different temperature 


boundary conditions. In all cases, it can be seen that 
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increasing PE, ie increasing dissipation, causes a reduction 


ripe nemload=Carrying Capacity of a, bearing. 
CONVECTION EFFECTS - NO INERTIA 


Figs. 47 through 49 present velocity distributions 
fOr be 1-0 7,.5-0, and 20.0, respectively ~ Upon examineng 
theirove locity wprofiles, we find that there as not. much 
difference among them, so that here the mass flow won't 
be important in explaining the variations in the pressure 
profiles. However, after examining the temperature distri- 
butions of Figs. 50 through 52, we find that the difference 
in the profiles is now quite pronounced, with the temperature 
profiles becoming smaller as ae increases. from, 1.0) to. 20.0. 

These results agree with findings made by a few 
previous researchers. Rodkiewicz and Anwar [11] found that 
LOmscelugqiePrandcl, NO.wilubricant, neglecting the convective 
terms in the energy equation caused a considerable rise in 
the lubricant film temperature; and a similar finding was made 
by snyder. (10). ) Both, groups concludeds thats forgnigneerandc) 
No. flows, it was necessary to retain both of the convective 
terms in the energy equation 

The pressure gradient is directly proportional to the 
lubricant viscosity, which increases with decreasing 
temperature. We should therefore expect that the pressure 
distribution curves would become larger as Pe increases. This 


expectation is proved correct by the results plotted on the 
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Again, since the viscosity is constant all along.the 
slider, then the shear stress at-the slider is proportional 


to the valueeot--the velocity gradient, ao at the ysaid 


interface. Figs. 47 through 49 show nen the velocity 
gradients. at the slider are all negative, and that the 
magnitude of the gradient at any section increases as Pe 
increases; so that the effect of increasing Pe is to increase 
the drag on the slider. Fig. 54 plots the shear stress 
distributions for the values of Pe under consideration, and 
it is seen that the results agree quite well with what 
we. expected. 

EVOS. Do. onreugh) S/eplot curves of loadmecapacicy. 
versus Pe for constant PE, but: for different temperature 
boundary conditions. In all cases,p1it can be seen that 


increasing Pe, ie increasing convection, causes an increase 


in the load carrying capacity of a bearing. 
SAMPLE CALCULATIONS 


In order to indicate the. physical significance of 
the dimensionless results, a sample calculation for the 
case of the cooled slider is presented below. This is 
followed by selected numerical results in tabular form. 


The following relationships have been used. 


Wu_UL 
Wwo= —3— (4.1) 
h 
O 
Du. UL 
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Soe fe Ws p Uh, (203) 
oC OH 
Pe = —P — (4.4) 
UG 2 
ig U 
Bee Cay Aan?) 
Bar 


5 lbf - Aaya 


#3 Btu s-sn 6 


lbf - sar - °F 


oeu were assumed 
sec - °F - in 


Following are the three cases for the cooled slider 


considered: 


(7) Assuming _T. .-@200EF ,-Tee= 1LO0°F, Pe = 2070 and 
s p 


PE = 1.0, we obtain T = (Te + Bee = 150°F and find 


respective absolute viscosity 


slider, 


Si lbf - sec 


2 


u. = 0.3276 x 10. 
£ an 


Having these, we find from (4.5) the velocity of the 


namely, 


U = 8.64 x 10° in/sec 


Now from the above given value of the Peclect number 


2 
3 Bak 


L 


Pe = 20.0 = 8.37 x 10 


gr x gor.t = = 1 bas 
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the value of Ao can be computed 


Doe 82 23) ax rome et 


Undenecnese conditions, from (4.1), (4.2), and. (4-3). we 
obtain 


111.2 Ibf/in- 


= 
I 


Ss 
Il 


aoe) ety. n (4.6) 


ee eee ligx 10°. yen 


(or) Fommvlanly for Wet=-50°R, T. = 100°F, Pe =) 2020 
s p 


and PE = 1.0 we obtain the following: 


feo tet 


FH 
I 


0.25798+x-10.% JE - sec 


ig ane 
in 


“Ss 
II 


ale Ta ee 107 in/sec 
=P: 
Ao = 4-684 x-lo al @) 


Under these conditions (4.1), (4.2) and (4.3) now yield 


Wa 21a Ibe /in- 

b= 1,945 1bf/in (4.7) 
2 fen 

Y= 2.77 x 10% ~ =e 


(i111) Assuming T = 5O°r, re = 100°F, Pe = 20.0 


: in, and using relationship (4.4) 
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Pe =. 20.0 = 8.37 x 10 
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the value of U can be computed 
ed 2°43 
U = 8.64 x 10° in/sec 


TY = 75°F, from which we find the respective absolute 


viscosity, namely, 


4 l1bf - sec 


2 


U, = 0.258 x eQ7 
in 


We can now compute the value of PE by using (4.5) 


UU 
PE = 61.25 = ct dey 8) 


1g 


Under the above conditions, from the relationships (4.1), 


(4.2) and (4.3) we now obtain 


Pago Eine 


WwW = 

Dt = -9.42 lbt/in (4.8) 
-2 lbm 

ee oO Det) LO Mn ae PS 


The above. format has been used to obtain the 


following numerical results: 


(a) | Slider cooled 
(b) Pad heated 
(c) Increasing PE 
(a) | Increasing Pe 
Wiile ithe, infitence of the inertia terms; onjthe pload 


capacity of a slider bearing is also investigated. 
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(e} Influence of Inertia Terms 
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SUMMARY OF CONCLUSIONS 


The main conclusions of the analysis may now be 
listed as follows. 

(1) Decreasing the slider temperature to pad 
temperature ratio of a slider bearing causes an increase 
in both the dimensionless load capacity and the dimensionless 
drag, irrespective of whether the slider is cooled or the 
pad is heated. For the range 0.5 = wee = 2.0; the 
dimensionless values of the load capacity, drag, and mass 
flow depend strongly on the es ratio, but weakly on the 
actual values of the temperature boundary conditions. 

(2) Although cooling the slider causes an increase 
in the load a bearing can carry, heating the pad for the case 
considered causes the bearing to be able to support less 
load. This difference is readily explained by using the 
following reasoning. 

For a bearing operating at a fixed speed and having 
a constant film thickness, the load is proportional to the 
product of the dimensionless load capacity and the reference 
viscosity. When the slider is cooled, both the values of 
the dimensionless load capacity and the reference viscosity 
increase, thus causing an overall increase in the load which 
the bearing can carry. However, for the heated pad, although 
there is an increase in the dimensionless load capacity, the 
value of the reference viscosity decreases, which for the 
case under consideration causes a drop in the load which 


the bearing can support. For both cases, the rate of mass 
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flow through the bearing increases, while for the cooled 
Slider alone there is.an increase inthe drag. 

(3) Neglecting the inertia terms in the momentum 
equation results in a smaller value being obtained for the 
load which the bearing can carry, and this error increases 
as Re- increases. The largest errors made in neglecting the 
inertia terms occur for the larger me’ te ratios. . Fer 
TS = 200°F, a = 1OO°F, and Re = 0.5, the value of the error 
is approximately 5% of the value obtained when the inertia 
terms are included, and this error rises to approximately 
10% of the true load for Re- = 3130 < 

(4) When the inlet to outlet ratio, hee of a plane 
slider bearing is increased, the load which the bearing can 
carry rises to a maximum then falls gradually. The observed 
Maxima were in tke range 2.00 ;\< hy Beeeey s 

(5) Increasing PE, brought about by varying the speed 
of the slider causes a drop in the dimensionless values of 
the load capacity and drag. 

However, with our choice of reference temperature, 
and for the sample example presented, the load the bearing 
Can carry increases about three times as PE goes from 1.27 
to 20.33. This result can be explained by using the following 
reasoning. When a bearing operates at a varying speed and 
for a constant film thickness, the load it can carry is 
proportional to the product of the dimensionless load capacity 
and the slider velocity assuming that the reference viscosity 


remains constant. For the larger value of PE, although the 
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value of the dimensionless load capacity decreases, the 
higher velocity offsets this decrease thus resulting in a 
larger value for the load which the bearing is able to 
support. A similar argument can be used to explain the 
larger drag. The increased dissipation caused by higher 
velocities would give a lower reference viscosity than the 
one used in our calculations, so that the actual- values of 
the load and drag would be smaller than those now evaluated. 

(6) Increasing Pe at constant speed, brought about 
by an increase in the film thickness results in larger 
values for the dimensionless load capacity and drag. 
However, the actual load the bearing can carry decreases as 
Pe increases, and the drag behaves likewise. The reason for 
this drop in load capacity can be understood from the 
following argument. When a slider bearing operates at 
constant speed and for a varying film thickness hoe the load 
tical carry 18° properctional toercie™ proauct =o. tne dimens1on— 


less load capacity and = , assuming the reference viscosity 


fo) 
remains constant. For the larger value of Pe, although the 


’ : ds 
value of the dimensionless load capacity increases, re 


fo) 
decreases because of an increase in ho: For the sample 


example presented this drop in = offsets the rise in the 
ho 

dimensionless load capacity, resulting in a smaller value 

of load which the bearing can sustain. A similar reasoning 


can be used to explain the smaller value obtained for the 


drag. The reduced dissipation caused by the larger film 
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thickness would give a somewhat larger value for the 
reference Viscosity, so that the actual magnitudes of the 


load and drag would be larger than those now obtained. 
SUGGESTIONS FOR FUTURE RESEARCH 


The present theoretical work has clearly shown that 
it 1s possible to considerably increase the load a slider 
bearing can carry by maintaining a low slider temperature to 
pad temperature ratio. It is suggested to design and 
build an apparatus to verify and find the range of 
applicability of the theory. The main feature of the design 
should be that a high viscosity layer next to the slider. 

One of the assumptions made in Chapter II was that 
the effects of thermal and elastic distortions were to be 
neglected in the analysis. However, for extreme cases 
examined, the difference in the temperature boundary conditions 
were so large that this assumption may not be strictly 
correct. The net result of thermal distortion would be to 
cause changes in the clearance between the slider and the 
pad, which directly affects the lubricant film thickness, 
and thus would cause changes in the load capacity, drag, etc. 
It is proposed to find the conditions at which these effects 


can still be neglected. 
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APPENDIX A 
TRANSFORMATION OF THE REDUCED EQUATIONS 


In order to simplify computations, the governing 
Equations (2.15) through (2.17) are transformed by use of the 


stream function. 


Writing eu = os , pvV=T- os and aves et Ss possible 


for us to obtain expressions for wu and v in terms of the 


stream function. 
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where we is the dimensionless mass flow. 
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APPENDIX B 
ENERGY EQUATION IN FINITE DIFFERENCE FORM 


mniewenergy Equation (2.25) is "rewritten 2n finvce 
difference form, by replacing the derivatives of the 
temperature by their central finite difference approximations. 
Let A x and Ay represent the grid spacing in the x and y 
directions as shown in Fig. 3. 
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APPENDIX C 
BEARING CHARACTERISTICS 


By substituting dimensionless variables for the 
dimensional ones, it is possible to obtain dimensionless 


equations for all the characterstics of the bearing operation. 
Dimensionless Mass Flow 
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APPENDIX F 
THE COMPUTER PROGRAMS 
Written in FORTRAN IV for an IBM 360/67 


Computer Using the MTS System 
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213 CONTINUE 
212 CONTINUE 
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GO TO 112 
i EVALUATION OF TEMP DIST 
ee = CSS a eee 


MM=MM+1 
25 NN=21 
ep RUN SS Se 
ITER=ITER+1 
DO 46 I=1,M1 
$$ _——————$.6—4° 7} 3=31 > 
BF (I,J) =G#HF (1) *AFM*SHY (1, J) /DX 
BE(T;3)= =BF (I, J) -2.D0*G 


raprayees 6/2. DO-HF (I) *APM*SHX (LL, 3) / (4. DO*DY) ) 


TP (3380: 2) W (1,3) =4B8F (1, 3) 

PESO SIE EI sti htt tt 
IF (J.#0.2)G0 TO 47 
@ (1,0) =4+BF(1I,3)—-ALF (I,J-1) *A(I,3) 


; byvj at 2 W {ay 
47 CONTINUE 
46 CONTINUE 
Oe ore rare eee a aaa ee 
DO 77 J=2,N1 
IF (J.EQ. NT) CNST=THETA(1+1, RES sad J)) 


IF (J. NE.2.AND.J.NE. N1)CNST=0, DO 
IF(NITER.LEQ. 1)B({1I,3)=THETA (1, 3-1) *({-A(I,J)) +THETA(I,J 


*THETA (I,J+1) * (-AH (I,J) ) #CEDRE* ( (SYY (I,J) #+SYY (I+1,J)) ¥*A 
*  EM/HP(I)) **2 
THEENS#O.1)909090 0 
IF (NITER.EO 1 TO $8 
LF (T (I,J) .GT. D2) VII=AB(I,J) x 
' iittala CEC Sire RT AE a 

IF(T (I,J). LB. 2. 3D2) VId=¥ (1,3) 
IF (T(I+1,d) .LB.2,3D2) VIL1=V (I+1,J) 
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ee sp peated a eo NOR ese e a aes ea nee ee 
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SYY (1, J) =H (1) **3¥*DPDX (I) *FY (I,J) /APM+DSUM (I,3)+C1(Z) *F 
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215 CONTINUE 
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214 CONTINUE 
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DESUM (I) =DESUM (1-2) +DX* (D(1) +4. DO*D (1-1) +D (I-2)) /3.D0 
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IF (ITER. EQ. 1.AND.IN.EQ. 0) MN=1 
IF (ITER. EQ. 1.ANDLIN.EQ.0)ITER=0 


IF (IN. £Q.1.AND. ITER. EQ. 2) MN=1 
IF (INJEQ.1.AND.ITER.GE.2)GO TO 269 
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IF (1. EQ. 2) EPVAL=DPVAL 
IF (1.EQ.2) GO TO 269 
eens eT ete a 
269 PB(I)=P (1) 
61 CONTINUE 
Sree 
IF(IN.EQ.0)GO TO 291 
IF (IN. EQ. 1.AND.ITER.GE.2)GO TO 291 
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TF (MN.EQ.0)GO TO 267 
‘i CLF PA HOKU CE Mes MOM At) 17 t= DOF Ty— —— 
UY (1, 1)= (4. DO*D (4, 2) -3,00"U (1, 1)-UX1,'3)9 (2 DOeDY) 
267 DO 108 J=1,N 
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IF (MN.EQ.0)GO TO 10 
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IF (MN.EQs0)GO TO 272 
IF (MNJEQ.1)GO0 TO 12 
SPEIRS ADS IED CE Ht oie ae 
IF (EPVAL.LT. EPRES) IN=0 
IF (EPVAL.LT.EPRES)GO TO 272 
EF 
e EVALUATICN OF AVERAGE GRAD'S FOR POINTS IN FILM 
272 werevieniws=+201 
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174 CONTINUE 
178 CONTINUE 
TEE of. 3 i Vere ae eee 
IF (EOVAL. LT. ETHET.AND.MM.EQ.1)G60 TO 12 
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139 FORMAT(8X,'RESULTS CONSIDERING VARYING VISCOSITY") 
TGTSty WRITE (6; 118) MITER; ITER, LOVAL — 

118 FORMAT(12X,'NITER EQUALS!,2X%,13,2X,*ITER EQUALS", 2x, 
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IF (EQVAL.LT. ETHEL) WRITE (6,119) ETHET 

119 FORMAT(40X,"EQVAL LESS THAN' , 3X,D13.6) 
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297 CONTINUE 
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DO 296 I=1,M 
WRITE (6, 301) (T(t (Oy , J=T, 42) 


296 CONTINUE 
301 FORMAT (2X, 12(F7.2,1X)/) 
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DEIDYP (I) =-H (1) *QYP (1) 
pI ters ts 
FDY (1) =AFM*¥SYY (1,1) / (5 (1) **2¥*F (1, 1)) 
13 CONTINUE 
aes so A oe be ie a nn eee SS 
195 FORMAT(30X,'*IEMP GRADIENTS AT PAD?) 
WRITE (6,197) (DTDYP (I) ,I=1,8) 
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PRINT 196 
196 FORMAT (30X,'TEMP GRADIENTS AT SLIDER‘) 
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PRINT 143 
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744 FORMAT (30X,'HEAT FLUX AT SLIDER") 
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124 ae YEES PETS FOR me cAe Fen DISTRIBUTION’) 
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WRITE (6, 33) (s(Z lat Be ea 
WRITE(6,161) (S(1, a3) con teee 


e EM? a3; 
161 FOR MAT (2X, 9{(F9. 4 1X) /) 
5071 CONTINUE 
BRENT —344+——_——— Sic cpoamale ee ee eee Daigle ami 
ZITOCFORBAT(S0X,*RESDLTS FOR YVELOCIT Dist*} 
BG3S02 DI= 1, 2s 
TE {(b;, 10) gS Ty 1. 
WRITE(6 7761) (U (1 7d) 55513727) 
502 CONTINUE 
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145 FORMAT (30X,'SHFAR STRESS AT SLIDER’) 
WRITE(6, 197) (EDY (I) p T= 1,48) 
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FSUM=0.D0 
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CPSUM=0.50 
OSSUM=0.130 
Dowd ~29511,2 
PSUM=PSUM+ (P (1) #4. DO*P (1-1) +P (1-2) ) *DX/3.D0 
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14 CONTINUE 
CSPER=DABS (QSSUM) / (DABS (QSSUM) +DABS (QPSUM) ) 
—5—-P RTT 0 ees 
140 FORMAT (30X,'FINAL RESULTS FOR PRESS DIST) 
adc 105) (P{I) ,1=1,6) 


RITE(6, 151) (P (1) ,I=7,11) 
151 Seda sa aes 5 (D13.6,2X)/) 
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PRINT 141 
141 FORMAT(2X,'LOAD CAPACITY!,10X,"DRAG',10X,'MASS FLOW! 
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PTHRU SLIDER") 
WRITE (6,142) PSUM,FSUM,AFM, QSPER 
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WRITE (6,166) OSSUM, OPSUM, FSUM 
5 FORMAT SE Be 6 Ab SS ve ama ee 
pO.199 J=1,N 


e: EVALUATICN OF HT CONVECTED INTO BRG 


ee * : y . 7 
e EVALUATICN OF HI CONVECTED OUT OF BRG 
TSEM (J) =APM*THETA (M,J) *SY (M, ,J) 
—+38-forttye2 ee 
OISM=0. D0 
QESM=0.D0 
ere ment heC Ts sek, eee 
OISH=QISM+DY* (TSIM (J) +4. D0* TSIM (J-1) +TSIM(J-2)) /3.D0 
QESM=QESM#DY* (TSEM (J) +4. DO*TSEM (J-1) +TSEM (J-2)) /3. DO 
UR SesT ie arenas aed a 
WRITE (6,201) QISM,AFM 
201 FORMAT(4X,"HEAT CCNVECTED INTO BRG EQU?,2X,F10.5,? 
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@RITE(6, 202) OESM, AFM 
202 FORMAT(4X,'HT CCNVECTED OUT OF BRG EQU',2X,F10.5, 
See 
c EVALUATION OF AVERAGE TEMP IN BEARING 
: DO 799 I=1,M 
——— he aNaGy 0p —guwrt SS 
UTSM=0.DO 
TSM=0.D0 
EEOC 8 ee aa Se 
USM= OSM+DY* (U (I,J) +4. DO*U(I,J-1) +U(I,J- 2)) /3. DO 
UTSM=UTSM+DY* (U (I,J) *THETA (I,J) +4.D0* (U (I,J-1) ¥THETA (I 
—— iUEe CGR Ree h) oo So ee ee SS 
*4+0 (I,J-2) *THETA (1,9-2)) /3.D0 
ISM=TSN+DY* (THETA (I,J) +4.DO*THETA (I,J-1) #THETA (I,J-2)) 
ee —— ; 
798 CONTINUE 
TB (I) =TAV*UTSM/USM 


799 CONTINUE 
WRITE (6,105) (TB(I) ,I=1,6) 
$$ ge tre tet tytn (ye ete 
WRITE (6, 105) (TA(I) ,I=1,6) 
WRITE (6,151) (TA(1I) ,I=7,11) 
TASM=0.DO 
POWFOTLE=3, NY2 
- Trea ST et st 0 OO 
TASM= TASM4+DX* (TA (1) 44. DO*TA (I-1) +#TA {I-2)) /3.D0 
797 CONTINUE 
+ h$$#+6,764) 4484, 4S 
784 FORMAT(' ',/,/,30X,*BULK TEMP EQUALS',2X,E£13.5,2x, 
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* EQUALS',2X,E13.5) 
Ft £0s 1) 6010-6 3 
NA=NA+1 


TP (NA. EQ. IVPRINT277 


a. 


; : - 
IF (NA. EQ.1. AND. IN.EQ.0)G0 TO 275 
ANS (I1,1)=R 

es oS 
ANS (11,3) =T! 
ANS {11,4)=ISTP 
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ANS (11,6) =PEC 
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ANS (11, 10)=CFDRE 
HO il, = | 
ANS {11, 12) =FSUM 
ANS (211,13) =AFM 
oes See eae 
ANS (11,15)=0SS 
ANS({11, a Behe 
p = SH | 
ANS (11, 18) =BT 
ANS (11,19) =CNTHT 
amc <TR GR PTR A ca ama a ge a ke 
ANS (11,21) =TES™ 
ANS{11,22)=TASM 
ee a Wears a ee ee pe ee ee ee ee eee 
PRINT 698 
698 FORMAT(*1!,/////,30X,*RESULTS OF NUMERICAL ANALYSIS 
ee ee PE amen eats a 
DO 697 11=1,NU8 
WRITE (6,700) (ANS (11,3) ,3=1, 7) 
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WRITE (6, 899) {ANS {21.3} ,3-8 715) 
899 FORMAT {(* ',3%,7{2173.5,2%))} 
" ee en eee — 
Hu FORMAT (* 1 3% ,8t213.5,22)) 
691 CONTINUE 
eee Tee ese 2 ee ee ee ee ee a a arr ae 
END 
SUBRCUTINE INV 
OT TOT TALE NED gg OL 
DIMENSION 1T(67),W(67),ARG(10) ,VAL(10) 
REWIND 2 
em SB de eee ee ee pee ee ee a Cr ee 
RETURN 
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ENTRY VIS(TEMP,V) 
RET Bao pis tO is hee, ko; OY 
iF (TEMP. LE. 2. 3D2) ERR=0.0168E-6 
IF (TEMP.LE.1.1D2) ERR=0. 148-6 
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CALL DALI(TEME,ARG,VAL,V,10,ERK,IER) 
RETURN 
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